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Abstract 



Using the supersymmetry approach, we study spectral statistical proper- 
ties of a two-dimensional quantum particle subject to a non-uniform magnetic 
field. We focus mainly on the problem of regularisation of the field theory. 
Our analysis begins with an investigation of the spectral properties of the 
purely classical evolution operator. We show that, although the kinetic equa- 
tion is formally time-reversible, density relaxation is controlled by irreversible 
classical dynamics. In the case of a weak magnetic field, the effective kinetic 
operator corresponds to diffusion in the angle space, the diffusion constant 
being determined by the spectral resolution of the inhomogeneous magnetic 
field. Applying these results to the quantum problem, we demonstrate that 
the low-lying modes of the field theory are related to the eigenmodes of the 
irreversible classical dynamics, and the higher modes are separated from the 
zero mode by a gap associated with the lowest density relaxation rate. As a 
consequence, we find that the long-time properties of the system are charac- 
terised by universal Wigner-Dyson statistics. For a weak magnetic field, we 
obtain a description in terms of the quasi one-dimensional non-linear a-model. 

I. INTRODUCTION 

Over recent years, the phenomenon of quantum chaos has been the subject of intense 
theoretical and experimental investigations |T[. Of these, the most recent has been the de- 
velopment of a statistical field theory in which the quantum properties of classically chaotic 
systems are expressed through the modes of density relaxation specified by the classical 
evolution operator 0. In the present paper, we will apply this field theoretic procedure 
to investigate the quantum dynamics of a two-dimensional particle subject to a random 
magnetic field (RMF). Our choice is motivated by two factors. Firstly, this problem rep- 
resents one of the simplest examples of a chaotic system. Secondly, it is of great practical 
interest, and as such has already received a lot of attention in the literature P-jlTf. To put 



the present investigation in perspective, we will first describe the basis of the general field 
theoretic construction. Then, we will review the main results on the RMF problem, and 
outline the strategy for the rest of the paper. 

The field theoretic approach to general chaotic quantum structures has been motivated 
by the success enjoyed by the statistical field theory of weakly disordered metallic conduc- 
tors. More specifically, applied to a model in which time-reversal symmetry is broken, an 
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average over realisations of a Gaussian (^-correlated random impurity potential shows that 
the two-particle spectral properties of a weakly disordered Hamiltonian are described by a 



super symmetric non-linear a-model with the effective action [JTS 



S - 

4 



J drstr [hD (VQ) 2 + 2is + Aa* R Q] , (1) 



where v = 1/VA denotes the density of states per unit volume V, D defines the clas- 
sical diffusion constant, and sA represents the symmetry breaking energy source with 
(Tg 15 = diag(l, — 1) A r- Here Q(r) represents a 4 x 4 supermatrix field obeying the non- 
linear constraint Q 2 {r) = 1. (An extensive review of this standard formalism can be found 
for example in Ref . ]18| . ) 

The connection between the supersymmeric field theory and physical coherence effects 
in disordered conductors can be understood by analysing characteristic scales of the theory. 
The effective action ([!]) identifies two time scales, the diffusion time to = L 2 /D and the 
Heisenberg time tn = ti/A. Their ratio defines the dimensionless conductance which in 
a good metal is large: g = tuj^D 3> 1- For energies sA <C h/tp, where quasi-classical 
dynamics is ergodic, the effective action is dominated by the zero spatial mode, Q(r) = Qo, 
independent of r 

S[Qo] = str K R Q ] • (2) 

In this limit, spectral properties are universal and, in this case, coincide with those of ran- 
dom matrix ensembles of unitary symmetry. Conversely, on energy scales sA ^> h/tp, the 
effective action is dominated by the global saddle-point Q sp = a^ R . An expansion in terms 
of the generators of the coset identifies the low-lying modes of density relaxation as diffu- 
sion modes. Interaction of these modes induces the well-known quantum weak localisation 
corrections. 

Surprisingly, an analogous theory can be defined for systems which are non-integrable 
but not stochastic 0. The key to the construction of an effective field theory of individual 
chaotic structures relies on the recognition that the spectrum itself provides a statistical 
ensemble By averaging over a wide interval of energy iVA it is possible to develop a 
"ballistic" non-linear a-model within the same general framework. In this case the effective 
action takes the general form 

S[Q] = -i| J dxystr [s + af Q + 2ihT- l CT] . (3) 

where £ = {H, } denotes the classical evolution operator, and the supermatrix field Q = 
T~ l <r^T depends on the 2d — 1 phase space coordinates xii = (r, p)2d-i parametrising 
the constant energy shell H(r, p) = E. (Here the coordinates are normalised such that 
/ cfxn = 1.) In the presence of a magnetic field, Q-matrices in (0) correspond to the unitary 
ensemble and are subject to additional symmetry constraints - see Ref. |20|] for details. 



While the zero-mode contribution to the action (H) reproduces random matrix or Wigner- 



Dyson statistics [21]], higher mode fluctuations establish non- universal quantum corrections. 



The utility of the field-theoretical approach beyond the universal regime was demonstrated 
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in papers by Blanter et al. p2| |, and Tripathi and Khmel'nitskii [^3| where the effective 



action (||) was applied to the study of a quantum billiard with diffusive surface scattering, 
a system that exemplifies a ballistic system in the regime of strong chaos. However, the 
cr-model approach has so far failed to provide explicit results for truly ballistic systems. 
This can, in part, be attributed to the fact that the kinetic operator £, which enters the 
effective action, is antihermitian for any ballistic system, a signature of the reversible nature 
of classical dynamics. As a consequence, its eigenvalues lie on the imaginary axis, so that 
the time evolution of a distribution in the phase space does not exhibit relaxation into the 
uniform ergodic state. Moreover, for chaotic systems, £ is ill-defined in the following sense: 
Classical dynamics involves stretching along the unstable manifold and contraction along 
the stable one. Thus, any non-uniform initial distribution will evolve into a highly singular 
function. In terms of the field theory, it means that the functional integral suffers from 
ultraviolet divergencies. In order to understand this, let us notice that the action (|3|) is only 
sensitive to the variations of the Q-matrix along the classical trajectories. Therefore, nothing 
prevents the Q-matrix from fluctuating in the transverse direction. It is these short-scale 
fluctuations that ultimately lead to the divergence of the functional integral |19[ . 



Several regularization procedures have been proposed to circumvent this problem [24-27 



One of the most natural ways to regularize the functional integral is to introduce an addi- 
tional term into the effective action which suppresses the fluctuations of Q in the transverse 
direction. After performing the integration, one should take the regulator to zero. A surpris- 
ing feature of the chaotic dynamics is that the limits time-to-infinity and regulator-to-zero 
do not commute. In particular, the eigenvalues of the regularized kinetic operator in the 
limit regulator-to-zero remain complex, with finite real parts corresponding to relaxation 
rates into the equilibrium distribution. These complex eigenvalues, which are independent 
of the regularization procedure, reflect intrinsic irreversible properties of the classical chaotic 
dynamics, and are known as Ruelle resonances or the Perron- Frobenius spectrum [2^,29]. 

One should stress that the direct regularization method described above is of little prac- 
tical use. This has been the main reason why all the attempts to apply the ballistic cx-model 
to real systems have been unsuccessful. In the present paper we propose a new general ap- 
proach to the regularization problem, which allows one to construct an effective field theory 
for the low-lying part of the Perron-Frobenius spectrum, and is applicable for both classical 
and quantum systems. To describe the basis of this approach, let us recall that the kinetic 
operator is defined in some Hilbert space, elements of which are smooth functions of the 
phase space coordinates. As was mentioned earlier, the eigenf unctions of £ are highly sin- 
gular and lie outside the Hilbert space. However, by properly choosing a subspace of the full 
Hilbert space one can make £ a well-defined operator in the sense that its eigenfunctions 
belong to the same subspace. Obviously, this subspace should correspond to the physically 
relevant low-lying part of the Perron-Frobenius spectrum. As we will see later, the kinetic 
operator defined in such a way is irreversible. To summarize, instead of calculating the full 
spectrum of the regularized kinetic operator, one can instead obtain an effective operator 
that correctly describes the low-lying modes of the irreversible classical dynamics. 

In order to choose the subspace in which the effective kinetic operator will act, one should 
truncate the basis in the Hilbert space by eliminating degrees of freedom irrelevant for the 
long-time evolution. Operationally, this is best done in the field-theoretical formulation. To 
illustrate our approach, let us consider the following formalism, wherein the Green function 
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of the classical kinetic operator is represented as an integral over the superfield if) = (i/)b, ^f), 



g(u) = - — - — j = / Dtp* Dip ipB ® %p B exp 



(4) 



iuo — C 

We then employ an RG-type approach and decompose ip in the following way, 

^(x||) = *(x||)+ X (x||) ) (5) 

where ^(xy) denotes "slow" fields, and x( x ||) "fast" fields. Integrating over the "fast" fields 
X one obtains an expression for the effective Green function G in the \I/-subspace, 



j D^*Dty ty B ® V* B exp - /'t2x||* t (x||)(iw-^ eff )*(x||) , (6) 



containing a renormalized kinetic operator C e g. There is no universal recipe as to how one 
should identify the "fast" and "slow" fields. However, this can be done straightforwardly in 
the case of "weak non-integrability"[], if the kinetic operator can be represented in the form 

C = C-0 + £pert; (J) 

where Co describes an integrable classical system, and C peTt can be considered as a small 
perturbation. In this case one can, in principle, calculate the eigenfunctions and eigenfre- 
quencies of the unperturbed operator Cq. "Fast" fields are identified as eigenfunctions of Co 
corresponding to high frequencies, while "slow" fields correspond to the low-frequency part 
of the spectrum of C . 

The main purpose of the paper is to demonstrate how this general regularization pro- 
cedure works by applying it to the problem of a particle confined to two- dimensions and 
propagating in a random magnetic field. The latter has attracted great interest in recent 
years. Firstly, there exist a number of experimental realizations where a RMF is imposed on 
2D electron gas. These include cases where the RMF is imposed by randomly pinned flux 
vortices in a type-II superconducting gate 0, by grains of a type-I superconductor ||, or by 
a demagnetized permanent magnet placed in the vicinity of the electron gas ||. Secondly, 
models of this type have been proposed within the gauge theory of high-T c superconductiv- 
ity [§,0. Thirdly, the problem is thought to be relevant for the composite fermion theory of 
the quantum Hall effect near v = 1/2 HfJ. 

The classical description of the transport properties of two-dimensional charged parti- 
cles in a random magnetic field with long-range correlations can be given in terms of the 
distribution function /(r, p), which obeys the Boltzmann equation 

OT OT C Op T 



To eliminate possible confusion, we stress that the system is assumed to be fully ergodic through- 
out the paper; in other words, we are in the regime of hard chaos, and the KAM theorem is not 
applicable here. 



4 



where the right-hand side describes relaxation to the uniform distribution / = 8(v — Vf) 
due to impurity potential scattering. Assuming the magnetic field is weak, one can calculate 
the averaged Green function of Eq. (RJ) using a perturbation theory. Within the Born 



approximation the renormalized value of the transport scattering time r tr can be found [10 



In particular, in the absence of the potential scattering (r — ► oo), one obtains for the 



transport time [11 
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/ dr(B(0)B(r)), (9) 
Jo 



r tr Tripp Jo 

where the angle brackets denote the ensemble average. 

In the opposite case of a strong magnetic field, where the Born approximation is not ap- 
plicable, a number of other approaches have been developed for solving the Boltzmann equa- 
tion, including the "Eikonal" approximation | 12| , and the percolation theory [p[| . Another 



possible approach is to use known results from the time-independent diffusion-advection or 
Passive scalar equation 

= DV 2 n + u • Vn, (10) 

at 

where u is a random velocity field. Indeed, starting with Eq. (RJ) and assuming a rapid 
relaxation in the momentum space, one obtains an effective equation for the particle density 
n(r) = J dpf(r, p), which has the form of Eq. (|T0|). It can be interpreted as the advection 
of guiding centers (or "van Alfven drift" of cyclotron orbits) due to the random component 
of the magnetic field. The diffusion-advection problem has been extensively studied because 



of its importance in fluid dynamics and plasma physics ||13|| . Known results suggest a 
non-trivial scaling behaviour of the conductivity a xx , which has been confirmed in recent 
experiments [Q . It is also interesting to note that an analogy exists between the diffusion- 
advection problem and the localization problem for a quantum particle in a random vector 
potential [p~5|[ . 

The quantum transport properties of a particle subject to a random magnetic field have 
been studied by Aronov, Mirlin, and Wolfle Jl6| . Applying a Gaussian 5-correlated distribu- 
tion for the magnetic field, spectral and transport properties are shown to be described by a 
diffusive supersymmetric non-linear cx-model of unitary symmetry. Although the relaxation 
rate which enters the average single-particle Green function is divergent if calculated pertur- 
batively, the transport relaxation time, which is a characteristic of the average two-particle 
Green function, turns out to be finite. The question of finding a physically meaningful 
definition of the single-particle relaxation time has been discussed by Altshuler et al. [|Hj[. 

In contrast to the above work by Aronov et al, we want to focus on the case of a 
weak slowly-varying magnetic field. It is convenient for our purposes to divide the rest of 
the paper into two sections. In the first section, we deal with the classical problem. The 
procedure of separating slow and fast fields, outlined above, is carried out in detail, and 
the effective classical kinetic operator is evaluated. The latter is shown to contain a term 
describing diffusion in the angle space. We calculate the corresponding "diffusion constant" 
(transport time), which is determined by the power spectrum of the magnetic field. The 
angle diffusion term is responsible for the relaxation of the momentum- dependent degrees 
of freedom, leading to the conventional diffusive dynamics at large scales. 
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In the second section, we consider the quantum problem, and apply the quasi-classical 
field theoretic description (^j. We show that, cast in the form of a field theory, the classical 
and quantum problems are closely related. One can therefore separate slow and fast degrees 
of freedom in Eq. (||) in almost the same way as one does for the classical problem. The 
only complication is the somewhat more complex structure of the integration manifold. 
Taking this into account, and using results from the first section, we arrive at the following 
conclusions: At short scales (small system size), spectral properties are described by a quasi- 
one-dimensional a-model. At intermediate scales, we obtain a ballistic a-model, where the 
traditional collision integral is replaced with diffusion in the angle space. At very large 
scales, we recover the conventional diffusive a-model. In this way, we demonstrate that a 
description in terms of the diffusive cr-model holds not only for systems with 5-correlated 



magnetic fields dominated by quantum scattering JRl , but also for semi-classical motion in a 



long-range magnetic field. In other words, although the semi-classical description is formally 
valid, quantum corrections are essential, and the transport at large scales is dominated by 
localization effects. Finally, it is important to emphasize that, in contrast to the existing 
literature on this problem, we do not make use of the notion of an ensemble average, so 
that our results are valid for an individual system with a given configuration of the magnetic 
field. 



II. CLASSICAL DYNAMICS 

As mentioned above, as a precursor to the analysis of the quantum problem, we begin 
by considering the dynamics of a classical particle in two-dimensions subject to an inho- 
mogeneous perpendicular magnetic field. Indeed, provided that h/mv <C a, where m is the 
mass of the particle, v is its velocity, and a is the typical length scale at which the magnetic 
field fluctuates, the motion of the quantum particle can be described semi-classically. Since 
energy is conserved, |v| = const., the motion is defined by the set of phase space variables 
(x, y, 0), which obey the following equations 

x = vcos<p, y = vsm(p, <j>=Q(x,y), (11) 

where Q = eB/mc represents the cyclotron frequency. We will assume that Q(x,y) is a 
random function of the spatial coordinates with zero average. We will further assume that 
the magnetic field is weak: Via <C v. 

As the particle moves through areas of different magnetic field, its trajectory is slightly 
deflected to one side or another depending on the sign of the magnetic field. From the point 
of view of the particle, the magnetic field Q(x(t),y(t)) is a random function of time. As 
such, one can obtain a good qualitative description of the particle's motion by taking this 
function to be 5-correlated Gaussian noise: 

<p = m, 

(Vi(t)Q(t')) = —5{t-t'). (12) 
From dimensional considerations it is clear that l/r tr ~ Q 2 a/v. From Eq. (|l2"D one obtains 
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2 (t)\ = f f dt'dt" (n(t')n(t")) = —t. (13) 

' JO JO T fr 



Such time dependence of (0 2 ) corresponds to diffusion in the angle space, with the diffusion 
coefficient l/r tr . If t ^ a/v, the motion is ballistic, and our description fails. At times 
a/v <C t <C r tr (the existence of this intermediate regime is ensured by the condition of 
weakness of the magnetic field fla <C v) one has ((f) 2 ) <C 1, so that the particle's trajectory in 
real space remains an approximately straight line. Finally, at times t > r tr , angle diffusion 
results in the uniform distribution of angles. Physically this means that the direction of 
movement becomes random, which leads to diffusion in real space (diffusive regime). Indeed, 
one has 

(x 2 (t)) = v 2 J* J* dt'dt" (cos cj)(t') cos 0(O) ~ 2D*, (14) 

where 

O = ^ (15) 

represents the diffusion coefficient in real space. 

This description offers a clear physical picture of what is going on in the system. How- 
ever, it has some drawbacks: Firstly, since (^-correlation contradicts the existence of a finite 
correlation length a of the magnetic field, it is not self-consistent. Secondly, a and fl them- 
selves are not well-defined quantities. Afterall, this is just a phenomenological description. 
We would like to obtain a quantitative description that would relate the diffusion coefficient 
to specific spectral properties of the random function Q(x, y). In order to do this we will 
make use of an alternative representation of the problem in the form of the Boltzmann 
equation. This representation has the advantage of being linear, allowing us to apply a 
field-theoretical approach. 

We start by recalling some basic facts about the Boltzmann equation: Any classical 
system obeying Hamilton's equations of motion can be alternatively described in terms of 
the distribution function /(x), which obeys the following partial differential equation 

^ + Cf = 0, (16) 

where the kinetic operator £ describes an incompressible flow in the phase space, 

d 

C = ±-—, <9-x = 0. (17) 
ax 

The incompressibility of the flow has a number of important implications, including the 
conservation of particle number, and the antihermiticity of C. In particular, for two arbitrary 
states |^) and \x), the following relation holds: 

(#C| X >* = -(xW). (18) 

Consequently, the eigenvalues of C are purely imaginary numbers. In the introduction we 
discussed the reasons why the spectrum of the kinetic operator for a chaotic system is not 
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well-defined, and came to the conclusion that the kinetic operator should be regularized. In 
the simplest case of a free particle, this can be done by introducing a term which accounts 
for the relaxation of the fast modes into the right-hand side of Eq. (|T6|): 



% + v-^f = ud 2 f, u>0. (19) 
at or 

Substituting / = f^exp^ik ■ r — iuit) into Eq. fll9D and taking the limit u — ■> 0, one finds 
that the eigenvalues iuQt) of £ acquire an infinitesimal positive real part: 

w(k)=v-k-i0. (20) 

As suggested in the introduction, the eigenvalues of a generic Perron-Frobenius operator 
have finite positive real parts. 

For the problem at hand, it is convenient to introduce the parametrization 

v = vn((j)), n (0) = (cos 0, sin 0). (21) 

The kinetic energy of a particle in the magnetic field is conserved, so that the phase space 
x = (r, p) is effectively reduced to the energy shell xy = (r, <fi), and the kinetic operator 
takes the form 

£ = ™(0).|- + fi(r)A (22) 

Since £ is a linear operator, Eq. ( |T^| ) can be solved. The general solution can be written 
in the form 

f{x,t) = U(t)f{x,0), (23) 

where U(t) = exp(— Ct) represents the classical evolution operator. It is often more conve- 
nient to work with its Fourier transform 



g(u) = - / dte**U(t) = j. (24) 

If we had an ensemble of different configurations of the magnetic field, information about 
transport properties could be extracted from the ensemble averaged distribution function 
(/(x, £)). As follows from Eqs. (^) and (H), it would be sufficient to calculate the ensemble 
averaged evolution operator (U(t)) or the Green function (g(u)). One could then introduce 
the effective "ensemble averaged" kinetic operator £ e g in the following manner: 

(U(t))=exp(-C cS t). (25) 

Since C is antihermitian, U(t) is a unitary operator. The ensemble averaged unitary operator, 
however, is not necessarily unitary, which means that C c s is not antihermitian, so that the 
spectrum of £ e s has a real part. Completely reversible dynamics becomes irreversible after 
the averaging. 



8 



One can, in principle, derive the spectrum of relaxation modes by calculating (g(u>)). 
But, since we are dealing with an individual system, not with an ensemble of systems, we 
are bound to take a different route. Following the ideas explained in the introduction, we 
therefore represent the classical Green function (p4|) as a field integral, 

g(u) = J DrD^i> B ®r B exv{S[rMi (26) 



where 



S = J dvd(j)^ (iu - 1) iff (27) 



denotes the effective action, and ip(r, 0) = (ipB, 4>f) represents a two-component supervector 
field, with indices B and F corresponding to bosonic and fermionic components. To make 
the functional integral (|26|) formally convergent, one usually puts an infinitesimal regulator 
into the action fl2"Tp. This is essentially the same regulator as in Eq. (EDI), and it will play 
an important role in the following analysis. 



The first term of the operator (|22[ ) describes the evolution of a free particle. Note 
that any function \I/(0) which depends on but not on r is necessarily a zero-mode of 
£free = vn((f)) ■ d / dr (i.e. £f iee ^((j)) = 0). This implies that the angle distribution (in 0- 
space) remains constant in time. Switching on a weak magnetic field, the angle distribution 
starts changing slowly. The evolution in 0-space is described by some effective operator 
which, following on from our naive physical picture, must have a diffusive form with some 
diffusion coefficient l/r tr . As we will see later, l/r tr is determined by spectral properties of 
0(r). 

This observation suggests a separation into slow and fast degrees of freedom from which 
the fast degrees of freedom can be integrated away. Since the particle's velocity is high 
and the magnetic field is weak, the spatial coordinate r(t) of the particle changes rapidly, 
while (p(t) changes slowly. The same statements are true for the distributions in the r-space 
and 0-space respectively. In terms of the action, one could say that spatially-dependent 
modes are fast, and spatially-independent modes are slow. The corresponding timescales 
are t c = L/v and r ir , where L is the size of the system. One can always choose L such that 
r c ^ T tr, and at the same time L^> a, the consistency of the two inequalities being ensured 
by the weakness of the magnetic field. (This is true only if the correlation length a of the 
magnetic field is finite - see below.) In this case, one can separate the modes with k = and 
perform the integration over the rest of the modes. By doing so, one obtains an effective 
action describing the evolution in the 0-space. For larger systems, a spatial dependence 
becomes important. The latter can be accounted for by dividing the system into smaller 
sub-systems and repeating the above procedure for each sub-system. Finally, in the case of 
very large systems, r c ^> r tr , the effective action is dominated by spatial diffusion modes. 

Beginning with small systems, we adopt the representation 

^(r, 0) = 4= + 4ff E Xk(0)e lk ' r , 0(r) = ^= £ fi k e ik ' r , (28) 



where Xk(0) are small fluctuating fields, and V denotes the volume of the system. Substi- 
tuting Eq. (EH]) into p7|), we obtain 
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r { r 1 d^f 1 



(29) 



where we have kept terms of up to second order. Performing the Gaussian integration over 
Xk(0), we obtain 



9W 

where 



1 J D**D** B ®** B exp{-S f &[**,*]}, (30) 



denotes the effective action, and 



|2 



< 32 » 

(Note that consistency requires only the zeroth-order term to be kept in the pre-exponential.) 
Replacing the summation in Eq. fl32|) by integration, and introducing the infinitesimal reg- 
ulator according to Eq. (GO), one finds 



1 . f \n k \ 2 dk 2tt f°° ,^ , 



— / dk\n k \\ (33) 



r fr j v • k — iO v Jo 

assuming that the spectrum of O(r) is isotropic, i.e. depends only on |k|. The effective 
action ([H]) corresponds to the kinetic operator of the form 

1 d 2 

£ cS = (34) 

T tr 0(j) 2 

The distribution function in 0-space /(</>) = J drf(r,<p) obeys the diffusion equation 

df 1 a 2 / 



<9t r tr 



S2 ' 



(35) 



Note that the functional integral (|30|) is convergent, since the eigenvalues of £ e fr are non- 
negative. In other words, C e g plays the role of a finite regulator in the field-theoretic 
representation. 

In conclusion, having started with a completely reversible kinetic operator with an imag- 
inary spectrum, we ended up with an effective irreversible operator describing a relaxation 
into the uniform distribution. It is important to emphasize that we did not introduce any 
collision integral. On the formal level, irreversibility comes from the infinitesimal imaginary 
part of the spectrum (p0[). One can draw an analogy between this problem and the problem 
of Landau damping in a plasma ||30|| . For the latter, the infinitesimal imaginary part of the 



spectrum ([2(]) results in a finite damping rate of the collisionless plasma. One could say 
that the initial infinitesimal relaxation rate (which one has to introduce simply to choose 
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the direction of the time evolution) is made finite by the regularization procedure, which for 
the problem at hand consists of truncating the phase space by eliminating large-k modes. 
This is quite a general situation: Many systems which exhibit irreversible properties such 
as increase of entropy and relaxation to equilibrium, are actually described by reversible 
classical mechanics. This phenomenon can be illustrated by the simple example of a gas of 
classical particles. A proper description must be in terms of the many-particle distribution 
function /(xi, . . . , xjy). The corresponding kinetic operator £(xi, . . . , xjv) is, of course, fully 
reversible, that is having a purely imaginary spectrum. In order to obtain a description in 
terms of the two-particle distribution function /(x, x') one has to truncate the phase space. 
As a result, the effective operator £ eff (x, x') contains a collision integral, which describes 
a relaxation into the Maxwell distribution. As shown by this example, there is nothing 
surprising in the fact that classical dynamics have irreversible properties. 

For the diffusion constant ( p3|) to be finite, there has to be a finite correlation length of 
the magnetic field. For a 5-correlated magnetic field one would have \Qk\ 2 = const., and the 
integral ( j3~3|) for l/r tr would diverge at large k, making the separation of modes illegitimate. 
We assume that the case a = corresponds to a different physical regime which can not be 
described by a diffusive approximation. 

Equation ( pM| ) is only valid if the distribution in the real space becomes uniform at times 
much shorter than the diffusion time in the angle space. This is the case if the size of the 
system L is not too large. In the opposite case of a very large system, L/v ^> r tr , one has 
to follow the spirit of the renormalisation group approach and divide the large system into 
smaller sub-systems. In each sub-system, one can separate the modes with k = and repeat 
previous calculations, integrating out the higher modes. The size of a sub-system should 
be small enough, but still much larger than the correlation length a of the magnetic field. 
Since each sub-system can be considered as independent, one has to introduce a coordinate 
R that labels the sub-systems: 



V»(r, <P) = -7=¥(R, 0) + -iff E Xk(R, 0)e 4k ' r , (36) 
VV V V wn 



;r v ' T ' . /t7 

where V stands for the volume of a sub-system. This representation effectively separates 
the degrees of freedom with small and large momenta. Substituting ( |36|) into (|27|), keeping 
the terms of up to the second order, and using the fact that the main contribution to the 
magnetic field comes from large momenta, we obtain: 



S = f dRd<p \ iu&V - J2 



+ 1 + 1 

it* • n(0) XkXk + ^ffX^— - W ^Xk 



(37) 



+ 1 



Integration over Xk(R, 4>) an d Xk(^> 0) yields 

1 0#t0vj/ 

which in turn implies a kinetic operator of the form 



r { , 1 4. 1 

s '« = ™+ - 1¥» - sri ' (38) 
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with the same transport relaxation time fl33|). 

We are most interested in the spectrum of the effective kinetic operator, which, as ex- 
plained in the introduction, coincides with the low-lying part of the Perron-Frobenius spec- 
trum. Representing an eigenfunction in the form 

M>(R,0) = ^(0)e 4q - R , (40) 

we obtain Mathieu's equation, which is formally equivalent to the Schrodinger equation for 
a quantum rotator in a uniform electric field, 

ip" + ivq cos (f>ip = iuip. (41) 



Equation (4~I) can not be solved analytically. However, if one is interested in the low-lying 
excitations, it can be recast in the form suitable for a perturbative treatment: 

1 d 2 

(H + V)?p = iujip, H = — V = ivq cos (/>. (42) 

r tr a<p 

The unperturbed eigenfunctions and eigenvalues are 

^) = -L e ^, «,,£) = — m 2 , m = 0, ±1, ±2, . . . . (43) 

V27T T tr 

Non-zero transition matrix elements exist only between neighbouring states, 

ivq 

'm-l,m *m,m—l ^ ■ (^ V 

Therefore, the first order corrections are zero. Calculating the second-order correction to 
the eigenvalue of the ground state (m = 0), we find 

Mq) = u>o = 1^ — —(&) = D( i ■ ( 45 ) 

m (^0) IUq — luJ rn 



The low-lying modes of the effective operator ( |39|) are nothing but diffusion modes in real 
space. The real space diffusion coefficient D is related to the transport time r tr by the usual 
equation flT5|). 

So far, we have shown that the effective Perron-Frobenius operator for the problem at 
hand contains a term that describes relaxation in the angle space. The effect of this term 
at large scales is to make the dynamics diffusive. Under the hypothesis of ergodicity (self- 
averaging), our formula (|33|) for the inverse scattering time matches the result obtained 
in the Born approximation after ensemble-averaging, Eq. (S). The advantage of the field- 
theoretical approach, however, is that it can be readily generalized to the quantum case, 
which is the subject of the next section. 
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III. QUANTUM DYNAMICS 



While the short-time dynamics of a particle in a slowly-varying magnetic field can for the 
most part be treated classically, the long-time dynamics are influenced strongly by quantum 
coherence effects which have manifestations both in localisation and spectral properties. To 
study the influence of these coherence effects on the quantum dynamics we will employ the 
statistical field theory defined in the introduction and applied to the quantum Hamiltonian 

£=^(p-^a) 2 , (46) 

where B = d x A. 

The action, defined by Eq. @, retains all the advantages of a simple classical description 
in the language of the kinetic equation, at the same time accounting for quantum interference 
effects. However, to apply the quasi-classical theory, we must introduce an appropriate 
regularization procedure. In doing so, we will take a lesson from the approach employed in 
the previous section. 

Let us first establish a connection between the quantum and the classical field theories, 
by making an expansion of the action (|3]) around the saddle point Q = a^ R , which is valid 
at large frequencies uj = sA/h. (In doing so, we make use of symmetry properties of T's - 
see Ref. [0.) Applying the parametrization 

T^ + i(° B 8 ) m , (47) 

expanding to the second order in B and B, and using the antihermiticity of C, we recover 
the classical action 



S[B, B]=Tiv J drd(f)stT [B (iu -£)&]. (48) 



At smaller frequencies, higher-order terms in the expansion generate quantum corrections. 

Assuming that the time scale at which the motion in real space becomes ergodic is set by 
the ballistic transport time across the system r c = L/v, we conclude that it is the spatially- 
dependent degrees of freedom that should be integrated out, in the same way as for the 
classical problem. Quantum effects will come from the interaction between the low-lying 
relaxation modes of the effective kinetic operator, which is induced by the non-linearity of 
the manifold of Q-matrices. 

In order to separate slow and fast degrees of freedom in the action @, instead of the 
linear decomposition (|5|) one uses the following non-linear decomposition, 

T(x||)=f(xn)T(x||), (49) 

where T(xy) stands for fast fields, T(xy) - for slow fields. It ensures that the Q-matrix always 
stays on the manifold Q 2 = 1. As explained above, fast fields are spatially-dependent, while 
slow fields depend only on angle 0. Since T(x.») fluctuates only weakly around T(x.i\) = 1, 
we arrive at the representation 
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T(t, 



E 



5_ k (0) 
5 k (0) 



„ikr 



(50) 



where T(0) describes spatially-uniform modes, B k (<p) and B k {4>) are small fluctuating fields 
describing the modes with non-zero momenta. Supermatrices T(0) and T _1 (0) parametrize 
the supermatrix Q(4>), 



Q 



(51) 



so that Q(4>) obeys the non-linear constraint Q 2 
notation 



1L. It is also convenient to introduce the 



M = Td (P QT- 1 . 
Using the relation [Q, d ( f ) Q} + = 0, one can prove the identity 

2str(M AR M RA ) = str(^g) 2 . 



(52) 



(53) 



Now we are prepared to derive an effective field theory for the problem at hand. Substi- 
tuting the representation ([)(]) into the action (|3]) with the kinetic operator (p2|), and making 
an expansion in B and B, we obtain 



S 



ivk ■ n((P)B k B k + ^B k Q k M AR + l -Q k M RA B k 



(54) 



ihuj 



J d^stiafQ. 



Integration over the fast degrees of freedom yields the effective action 

1 



eff 



h 
'8A 



difistr 



Ttr 



(d^QY + 2iua^ R Q 



(55) 



which has the form of a one-dimensional diffusive a-model, the diffusion taking place in 
the angle space. The transport time is given by the same formula ( j33|) as for the classical 
problem. To arrive at Eq. (0), we have used the identity fl53"|). 

If ujT tr ^> 1, the main contribution to the functional integral comes from around the 
saddle point Q = a^ R . In the opposite case, ur tr <C 1, one must take into account the 
whole supermatrix manifold. Spectral correlation functions are then determined by univer- 
sal Wigner-Dyson statistics. Corrections to this universal behaviour can be calculated by 
performing an expansion in 1/g analogous to that employed by Kravtsov and Mirlin ||31|| . 
where 
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At, 



(56) 



tr 



plays the role of dimensionless conductance. 

At larger scales the low-lying modes of the Perron- Frobenius operator acquire a spatial 
dependence. This can be accounted for by using another representation, instead of (|50"D, 
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that separates degrees of freedom with small and large momenta (compare with the previous 
section). The effective action takes the form of the ballistic cr-model, 



S cS = — [ drdfatr hr-VfW^) ■ VT - —(dsQ) 2 + iujafQ\ . (57) 
4 J L ' 2r tr 

Note that a similar type of action was proposed by Muzykantskii and Khmelnitskii JJ], 
where a collision integral was introduced as a phenomenological term to describe quantum 
scattering by a 5-correlated impurity potential. The second term in the action fl5"T|) de- 
scribes the diffusion in the angle space and effectively plays the role of a collision integral 
in the sense that it is responsible for relaxation of the momentum-dependent degrees of 
freedom. However, its origin in our case is quite different, since it appears as a result of the 
regularization of the field theory. 

Finally, one can establish a relation between the ballistic cr-model (|57j ) and the con- 
ventional diffusive cr-model using the procedure similar to that of Ref. 0. Anticipating a 
rapid relaxation in the momentum space, and a slow variation of the spatial modes, one can 
integrate out the momentum- dependent degrees of freedom to obtain the action (TT]) with 
the diffusion constant (|I5|). Therefore, the transport properties of the system at very large 
scales are dominated by the localization effects. 

To summarise, in this paper we applied the recently-developed field theory of quantum 
chaos to the problem of a particle propagating in a non-uniform magnetic field. We proposed 
a new approach to the regularisation of the ballistic cr-model (||]), and obtained an effective 
description of spectral and transport properties of the quantum system in terms of the low- 
lying modes of the irreversible dynamics of its classical counterpart. We argued that this 
method is quite general and can be applied to a variety of chaotic systems. 
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15 



REFERENCES 



[1] M. C. Gutzwiller, Chaos in Classical and Quantum Mechanics, (Springer- Verlag, New- 
York, 1990) 

[2] B. A. Muzykantskii and D. E. Khmelnitskii, JETP Lett. 62, 76 (1995) 

[3] A. Geim, S. Bending, and I. Grigorieva, Phys. Rev. Lett. 69, 2252 (1992) 

[4] A. Smith et al, Phys. Rev. B 50, 14726 (1994) 

[5] F. B. Mancoff et al, Phys. Rev. B 51, 13269 (1995) 

[6] L. B. Ioffe and A. I. Larkin, Phys. Rev. B 39, 8988 (1989) 

[7] N. Nagaosa and P. A. Lee, Phys. Rev. Lett. 64, 2450 (1990) 

[8] B. I. Halperin, P. A. Lee, and N. Read, Phys. Rev. B 47, 7312 (1993) 

[9] V. Kalmeyer and S.-C. Zhang, Phys. Rev. B 46, 9889 (1992) 
[10] Per Hedegard and Anders Smith, Phys. Rev. B 51, 10869 (1995) 
[11] A. D. Mirlin, D. G. Polyakov, and P. Wolfle, Phys. Rev. Lett. 80, 2429 (1998) 
[12] D. V. Khveshchenko, Phys. Rev. Lett. 77, 1817 (1996 
[13] M. B. Isichenko, Rev. Mod. Phys. 64, 961 (1992) 

[14] L. Zielinski, K. Chaltikian, K. Birnbaum, C. M. Marcus, K. Campman, and A. C. 

Gossard, Europhys. Lett. 42, 73 (1998) 
[15] Jonathan Miller and Jane Wang, Phys. Rev. Lett. 76, 1461 (1996) 
[16] A. G. Aronov, A. D. Mirlin, and P. Wolfle, Phys. Rev. B 49, 16609 (1994) 
[17] E. Altshuler, A. G. Aronov, A. D. Mirlin, and P. Wolfle, Europhys. Lett. 29, 239 (1995) 
[18] K. B. Efetov, Adv. Phys. 32, 53 (1983) 

[19] A. V. Andreev, B. D. Simons, O. Agam, and B. L. Altshuler, Nucl. Phys. B 482, 536 
(1996) 

[20] B. D. Simons, O. Agam, and A. V. Andreev, J. Math. Phys. 38, 1982 (1997) 
[21] O. Bohigas, M. J. Giannoni, and C. Schmidt, Phys. Rev. Lett. 52, 1 (1984) 
[22] Ya. M. Blanter, A. D. Mirlin, and B. A. Muzykantskii, Phys. Rev. Lett. 80, 4161 (1998) 
[23] V. Tripathi and D. E. Khmelnitskii, Phys. Rev. B 58, 1122 (1998) 
[24] P. Gaspard, G. Nicolis, A. Provata, and S. Tasaki, Phys. Rev. E 51, 74 (1995) 
[25] D. Ruelle, Statistical Mechanics, Thermodynamic Formalism, (Addison- Wesley, Read- 
ing, MA, 1978) 
[26] G. Nicolis and C. Nicolis, Phys. Rev. A 38, 427 (1988) 
[27] H. H. Hasegawa and D. J. Driebe, Phys. Rev. E 50, 1781 (1994) 
[28] D. Ruelle, Phys. Rev. Lett. 56, 405 (1986) 
[29] M. Pollicot, Ann. Math. 131, 331 (1990) 

[30] E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics, (Oxford, 1981) 
[31] V. E. Kravtsov and A. D. Mirlin, JETP Lett. 60, 656 (1994) 



16 



